Loop Quantum Gravity and Ultra High Energy Cosmic Rays 
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There are two main sets of data for the observed spectrum of ultra high energy cosmic rays (those 
cosmic rays with energies greater than ~ 4 x 10 18 eV), the High Resolution Fly's Eye (HiRes) col- 
laboration group observations, which seem to be consistent with the predicted theoretical spectrum 
(and therefore with the theoretical limit known as the Greisen-Zatsepin-Kuzmin cutoff), and the ob- 
servations from the Akeno Giant Air Shower Array (AGASA) collaboration group, which reveal an 
abundant flux of incoming particles with energies above 1 x 10 20 eV violating the Greisen-Zatsepin- 
Kuzmin cutoff. As an explanation of this anomaly it has been suggested that quantum-gravitational 
effects may be playing a decisive role in the propagation of ultra high energy cosmic rays. In this 
article we take the loop quantum gravity approach. We shall provide some techniques to establish 
and analyze new constraints on the loop quantum gravity parameters arising from both sets of data, 
HiRes and AGASA . We shall also study their effects on the predicted spectrum for ultra high energy 
cosmic rays. As a result we will state the possibility of reconciling the AGASA observations. 



I. INTRODUCTION 

In this article we are concerned with the observation 
of ultra high energy cosmic rays (UHECR), i.e. those 
cosmic rays with energies greater than ~ 4 x 10 18 eV. 
Although not completely clear, it has been suggested that 
these high energy particles are possibly heavy nuclei 0,0 
(we will assume here that they are protons) and, by virtue 
of the isotropic distribution with which they arrive to us, 
that they originate in extragalactic sources. 

A detailed understanding of the origin and nature of 
UHECR is far from being achieved; the way in which 
the observed cosmic ray spectrum appears to us still is a 
mystery and a matter of great debate. A first subject of 
interest (faced with the lack of reasonable mechanisms) 
is how such energetic particles have been accelerated to 
energies well above 4 x 10 18 eV by its sources. A sec- 
ond subject of interest is the study of their propagation 
in open space through the cosmic microwave background 
radiation (CMBR), whose presence necessarily produce 
a friction on UHECR making them release energy in the 
form of secondary particles and affecting their possibility 
to reach great distances. A first estimation of the char- 
acteristic distance that UHECR can reach before losing 
most of their energy was simultaneously made in 1966 
by K. Greisen and G. T. Zatsepin & V. A. Kuzmin 
U , who showed that the observation of cosmic rays with 
energies greater than 4 x 10 19 eV should be greatly sup- 
pressed. This energy (4 x 10 19 eV) is usually referred 
as to the GZK cutoff energy. Similarly and a few years 
later, F.W. Stecker calculated the mean life time for 
protons as a function of their energy, giving a more accu- 
rate perspective of the energy behavior of the cutoff and 
showing that cosmic rays with energies above 1 x 10 20 eV 
should not travel more than ~ 100 Mpc. More detailed 
approaches to the GZK-cutoff feature have been made 
since these first estimations. For example V. Berezin- 
sky & S.I. Grigorieva H, V. Berezinsky et al. and 



S.T. Scully & F.W. Stecker ||| have made progress in the 
theoretical study of the spectrum J{E) (i.e. the flux of 
arriving particles as a function of the observed energy E) 
that UHECR should present. As a result, the GZK cutoff 
exists in the form of a suppression in the predicted flux 
of cosmic rays with energies above ~ 8 x 10 19 eV. 

At present there are two main different sets of data 
for the observed flux J(E) in its most energetic sector 
(E > 4 x 10 18 eV). On one hand we have the observations 
from the High Resolution Fly's Eye (HiRes) collaboration 
group , which seem to be consistent with the predicted 
theoretical spectrum and, therefore, with the presence 
of the GZK cutoff. Meanwhile, on the other hand, we 
have the observations from the Akeno Giant Air Shower 
Array (AGASA) collaboration group [lfj, which reveal 
an abundant flux of incoming cosmic rays with energies 
above 1 x 10 20 eV. The appearance of these high energy 
events is greatly opposed to the predicted GZK cutoff, 
and a great challenge that has motivated a vast amount 
of new ideas and mechanisms to explain this phenomenon 
El El El El El 111 E3, 113 • H the AG AS A observations 
are correct then, since there are no known active objects 
in our neighborhood (let us say within a radius R ~ 100 
Mpc) able to act as sources of such energetic particles 
and since their arrival is mostly isotropic (without any 
privileged local source), we are forced to conclude that 
these cosmic rays come from distances larger than 100 
Mpc. This is commonly referred as the Greisen-Zatsepin- 
Kuzmin (GZK) anomaly. 

One of the interesting notions emerging from the pos- 
sible existence of the GZK anomaly is that, since ul- 
tra high energy cosmic rays involve the highest energy 
events registered up to now, then a possible framework 
to understand and explain this phe nomena could be of 
a quantum-gravitational nature pi Ebl Ell El El Ej| . 
This possibility is indeed very exciting if we consider the 
present lack of empirical support for the different ap- 
proaches to the problem of gravity quantization. In the 
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context of the UHECR phenomena, all these different 
approaches motivated by different quantum gravity for- 
mulations, have usually converged on a common path to 
solve and explain the GZK anomaly: the introduction of 
effective models for the description of high energy parti- 
cle propagation. These effective models, pictures of the 
yet unknown full quantum gravity theory, offer the possi- 
bility to modify conventional physics through new terms 
in the equations of motion (now effective equations of 
motion), leading to the eventual breakup of fundamen- 
tal symmetries such as Lorentz invariance (expected to 
be preserved at the fundamental level). These Lorentz 
symmetry breaking mechanisms are usually referred as 
Lorentz invariance violations (LIV's) if the break intro- 
duce a privileged reference frame, or Lorentz invariance 
deformations (LID's), if such a reference frame is absent 
[25L l26l l27j . Its appearance on theoretical as well as phe- 
nomenological grounds (such as high energy astrophysical 
phenomena) has been widely studied, and offers a large 
and rich amount of new si gna tures that deserve attention 

m in us m in in in in m . 

To deepen the above ideas, we have adopted the loop 
quantum gravity (LQG) theory El EH 1 one of the pro- 
posed alternatives for the yet nonexistent theory of quan- 
tum gravity. It is possible to study LQG through effec- 
tive theories that take into consideration matter-gravity 
couplings. In this line, in the works of J. Alfaro et al. 
[3^. El EH- the effects of the loop structure of space 
at the Planck level are treated semiclassically through 
a coarse-grained approximation. An interesting feature 
of these methods is the explicit appearance of the Plank 
scale l p and the appearance of a new length scale £^> l p 
(called the "weave" scale), such that for distances d <C C 
the quantum loop structure of space is manifest, while 
for distances d ^ C the continuous flat geometry is re- 
gained. The presence of these two scales in the effective 
theories has the consequence of introducing LIV's to the 
dispersion relations E = E(p) for particles with energy 
E and momentum p. It can be shown that these LIV's 
can significantly modify the kinematical conditions for a 
reaction to take place. For instance, as shown in detail 
in El|, if the dispersion relation for a particle i is (from 
here on, h = c = 1) 

El = A\y\ + m\ (1) 

(where Ei, pi and nii are the respective energy, momen- 
tum and mass of the zth particle, and Ai is a LIV param- 
eter that can be interpreted as the maximum velocity of 
the ith particle) then the threshold condition for a re- 
action to take place can be substantially modified if the 
difference 5 A = A a — At is non zero (a and b are two par- 
ticles involved in the reaction leading to the mentioned 
threshold) 42] . An interesting consequence of the above 
situation — for the UHECR phenomenology — is that the 
kinematical conditions for a reaction between a primary 
cosmic ray and a CMBR photon can be modified, lead- 
ing to new effects and predictions such as an abundant 
flux of cosmic rays well beyond the GZK cutoff energy 



(explaining in this way the AGASA observations). 

The purpose of this paper is to provide some tech- 
niques to establish and analyze new constraints on the 
LQG parameters (or any other LIV parameters), that 
will confidently rise when the experimental situation is 
clarified in a reliably way up to a certain energy scale. In 
the present case, and for the practical purpose of this pa- 
per, we shall assume that such energy scale is currently 
4 x 10 19 eV. Also, we shall attempt to predict (under 
certain assumptions) a modified UHECR spectrum aris- 
ing from the LQG corrections to the conventional theory, 
and consistent with the AGASA observations (though 
we shall analyze both, HiRes and AGASA sets of data, 
throughout this paper we will be more concerned with 
the possibility that the AGASA results are the correct 
one). To accomplish these goals, we have organized this 
article as follows: In section[Hj "Ultra high energy cosmic 
rays" , we shall give a brief self contained derivation of the 
conventional spectrum and briefly analyze it jointly with 
HiRes and AGASA observations. In section UTTl "Loop 
quantum gravity" , we will present a short outline of loop 
quantum gravity and its effective description of fermion 
and electromagnetic fields (relevant for the description 
of UHECR propagation) . In section IIVI "Threshold con- 
ditions" , we will analyze the effects of LQG corrections 
on the threshold conditions for the main reactions in- 
volved in the UHECR phenomena to take place. In sec- 
tionlVl "Modified spectrum" , we will show how the mod- 
ified kinematics can be relevant to the theoretical spec- 
trum J(E) of cosmic rays (we will present the obtained 
modified spectrum). Section IVTl "Conclusions", will be 
reserved for some final remarks. 



II. ULTRA HIGH ENERGY COSMIC RAYS 

In this section we will review the main steps in the 
derivation of the UHECR spectrum. This presentation 
will be useful and relevant for the description of the kine- 
matical effects that LQG corrections can have on the 
predicted flux of cosmic rays. The following material 
is mainly contained in the works of F.W. Stecker 
Berezinsky et al. 7] and S.T. Scully & F.W. Stecker 
0. 

A. General Description 

Two simple and commonly used assumptions for the 
development of the cosmic ray spectrum are: 1) that the 
sources are uniformly distributed in the Universe, and 2) 
that the generation flux F(E g ) of emitted cosmic rays 
from the sources is correctly described by a power law 
behavior of the form F(E g ) cx E g la , where E g is the 
energy of the emitted particle and 7 g is the generation 
index. 

One of the main quantities in the calculation of the 
UHECR spectrum is the energy loss —E~ 1 dE/dt. This 
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quantity describes the rate at which a cosmic ray loses 
energy, and takes into consideration two chief contribu- 
tions: the energy loss due to the redshift attenuation and 
the energy loss due to collisions with the CMBR photons. 
This last contribution depends, at the same time, on the 
cross sections a and the inelasticities K of the interac- 
tions produced during the propagation of protons in the 
extragalactic medium, as well as on the CMBR spec- 
trum. The most important reactions taking place in the 
description of proton's propagation (and which produce 
the release of energy in the form of particles) are the pair 
creation 



P + 7 



P- 



and the photo-pion production 

p + 7 — > p + 7T. 



(2) 



(3) 



This last reaction happens through several channels (for 
example the baryonic A and N, and mesonic p and u) 
resonance channels, just to mention some of them) and 
is the main reason for the appearance of the GZK cutoff. 



B. Some Kinematics 

To study the interaction between protons and the 
CMBR, it is useful to distinguish between three reference 
systems; the laboratory system JC (which we identify with 
the Friedman Robertson Walker (FRW) co-moving refer- 
ence system), the center of mass (cm.) system JC*, and 
the system where the proton is at rest JC' '. In terms of 
these systems, the photon energy will be expressed as to 
in JC and as e in JC' . The relation between both quantities 
is simply 



= 7<j(l — /3cos#), 



(4) 



where 7 = E/m p is the Lorentz factor relating JC and 
JC' , E and m p are the energy and mass of the incident 

proton, P = y/l — 7~ 2 , and 6 is the angle between the 
momenta of the photon and the proton measured in the 
laboratory system JC. 

To determine the total energy E£ ot = E* + e* in the 
cm. system, it is enough to use the invariant energy 
squared s = E? ot - pl ot (where E tot = E + uj and p to t 
are the total energy and momentum in the laboratory 
system). In this way, we have 



E tot = s = m 2 



2m rj 



(5) 



As a consequence, the Lorentz factor j c which relates the 
JC reference system with the JC* system, is 



E 



E 



7c 



{ml + 2m p e) 1 / 2 ' 



(6) 



Let us consider the relevant case in which the reaction 
between the proton and the CMBR photon is of the type 



P + l 



(7) 



where a and b are two final particles of the collision. The 
final energies of these particles are easily determined by 
the conservation of energy-momentum. In the JC* system 
these are 



77** 

£j n h — 



1 



2^/~s 



(8) 



Transforming this quantity to the laboratory system, and 
averaging with respect to the angle between the direc- 
tions of the final momenta, it is possible to find that the 
final average energy of a (or 6) in the laboratory system 
is 



< K,b >- 




(9) 



The inelasticity K of the reaction is defined as the 
average fractional difference K = AE/E, where AE = 
E — Ef is the difference between the initial energy E and 
final energy Ef of the proton (in a single collision with the 
CMBR photons). For the particular case of the emission 
of an arbitrary particle a (that is to say p + 7 — * p + a) , 
expression © allows us to write 



1 / mi- 



(10) 



where K a is the inelasticity of the described process. This 
is one of the main quantities involved in the study of 
the UHECR spectrum, in particular, when the emitted 
particle a is a pion. 



C. Mean Life t(E) 

To derive the UHECR spectrum it is imperative to 
know the mean life t{E) of the cosmic ray (or proton) 
with energy E propagating in space, due to the attenu- 
ation of its energy by the interactions with the CMBR 
photons. The mean life t(E) is defined through the re- 
lation 



t(E)- 



1 dE 
E~dT 



(11) 



col 



where the label "col" refers to the fact that the energy 
loss is due to the collisions with the CMBR photons. To 
explicitly determine the form of t(E), let us express l|ll|l 
in terms of the microscopic collision-quantities: 



r(E)- 1 = 



AE 1 
IT A? 



(12) 



where AE is the difference between the initial and final 
energies of the proton before and after each collision, and 
At is the characteristic time between collisions. Introduc- 
ing the inelasticity through its definition K — AE/E, 
and expressing the characteristic time in terms of the 
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scattering cross section and density p of the target pho- 
tons, we can write then: 



r(EY 



Kap Urol, 



(13) 



where v rc \ is the relative velocity between the incident 
proton and the background. The above relation can be 
driven to a more accurate version if we consider that 
both v re \ and a are functions of the energy and direction 
of propagation of the CMBR photons relative to the in- 
cident proton. Considering these elements we are able to 
write 



dr(E) 1 = K a Vr C iT](u>)cLj dQ/Air, 



(14) 



where rj(uj)du> is the CMBR density of photons with ener- 
gies in the range [uj, Lu + dio], and dQ/Air = sm.6 d9d(f>/Air 
is the section of solid angle. With the above quantities, 
it is simple to rewrite v re \ through 



v Te i dQ /47r 



; de < 



Attj 2 ui 



(15) 



with e € [0, 2juj] and <fi 6 [0, 2-k\. Substituting equa- 
tion HHJ) in JT3J and using the fact that the CMBR 
density corresponds to a Planck distribution rj(ui)duj = 
L0 2 duj /7r 2 (e"/ fcT — 1), it is finally possible to show that 
the mean life t(E) can be written in the form 



r(E)- 1 = 



kT 



°°dea{e)K{e)e\n[l - e - c/27fcT ].(16) 



D. Energy Loss and Spectrum 

The energy loss suffered by a very energetic proton 
during its journey, from a distant source to our detec- 
tors, is not only produced by the collisions that it has 
with CMBR at a particular epoch. There will also be 
a decrease in its energy due to the redshift attenuation 
produced by the expansion of the Universe. At the same 
time, such expansion will affect the collision rate through 
the attenuation of the photon gas density, which can be 
understood as a cooling of the CMBR through the rela- 
tion T = (1 + z)Tq, where z is the redshift and To is the 
temperature of the background at the present time. To 
calculate the spectrum we need to consider the rate of 
energy loss during any epoch z of the Universe. 

For the present discussion, we shall assume that the 
Universe is well described by a matter dominated Fried- 
man Robertson Walker (FRW) space-time, and that the 
ratio of density fir, = p/p c (where p is the energy den- 
sity of the present Universe and p c is the critical energy 
density for the Universe to be flat) is such that Oq = 1- 
The above assumptions give rise to the following relation 
between the temporal coordinate t (proper time in the 
co-moving system) and the redshift z: 



dt = - 



dz 



H {l + z) 5 / 2 



(17) 



where Hq is the Hubble constant at present time. Since 
the momentum of a free particle in a FRW space behaves 
as p oc (1 + z), we will have, with the additional consid- 
eration p» m (where m is the particle mass), that the 
energy loss due to redshift is 



(18) 



On the other hand, the energy loss due to collisions 
with the CMBR will evolve as the background tempera- 
ture changes (recall T — (1 + z)Tq). This evolution can 
be parameterized through z and is given by 



1 dE 
~E~dl 



^37 ={l + zYr{[l + z]EY 



(19) 



col 



The total energy loss can be expressed as the addition 
of the former contributions (using z instead of t) 



1 dE 
~E~dz~ 



(1 + z)- 1 + H-\l + z)^ 2 t([1 + z]E)-\ (20) 



Equation (|2U[) can be numerically integrated to give the 
energy E g (E, z) of a proton generated by the source in a 
z epoch and that will be detected with a energy E here 
on Earth. Let us designate this solution by the formal 
expression 



Eg(E, 



X(E, z)E. 



(21) 



It is also possible to manipulate equation l|20l) to ob- 
tain an expression for the dilatation of the energy interval 
dEg/dE. To accomplish this it is necessary to integrate 
(|20|l with respect to z and then differentiate it with re- 
spect to E to obtain an integral equation for dE g /dE. 
The solution of such an equation is found to be 



dE g( z g 
dE 



x exp 



(1 + z g ) 



I' 1 



z) 



1/2 



db{E') 



uo 



dE' 



(22) 



where E' = (1 + z)X{E, z)E. 

The total flux dj(E) of emitted particles from a volume 
element dV — R 3 (z)r 2 drdH,, in the epoch z and coordi- 
nate r, measured from Earth at present with energy E 
is 



dj{E)dE 



F(E ,z)dE n(z)dV 
(1 + z)47ri?2 r 2 ' 



(23) 



where j(E) is the particle flux per energy, F(Eq, z)dEo 
the emitted particle flux within the range (Eq, Eo + dEo), 
and n(z) the density of sources in z. As previously men- 
tioned, it is convenient to study the emission flux with 
a power law spectrum of the type F(E) oc E~ lg . It can 
be shown that with such assumption, the relation be- 
tween the emission flux and the total luminosity L p of 
the source is F(E) = (j g — 2)L p E~ lg . To describe the 
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evolution of the sources we shall also use a power law 
behavior. This will be done through the relation 



L p (z)n(z) = (l + *)( 3 +™)L p (0)n(0), 
= (l + z)( 3 +™)£ , 



(24) 



in such a way that m = corresponds to the case in 
which sources do not evolve. If we consider that Rq — 
(1 + z)R(z) and R(z)dr = dt for flat spaces (and v ~ 1 
for very energetic particles), using (|17|l to express all in 
terms of z and, finally, integrating ill'-il) from z = to 
some z — z max for which sources are not relevant for the 
phenomena, it is possible to obtain 



m = (7.-2)^1^- 



,dEg{Zg) 



< y dz s (l + z s ) m - 5 / 2 A^(i?,z g )- 



(25) 



The above expression constitutes the spectrum of 
UHECR. It remains to fix (observationally) the volumet- 
ric luminosity Cq and the 7 g and m indices. 



A HiRes I 
■ HiRes II 




10 ,B 10 ia 10 2U 10 2i 

Energy E [eV] 

FIG. 1: UHECR spectrum and HiRes observations. The fig- 
ure shows the UHECR spectrum J(E) multiplied by E 3 , for 
uniform distributed sources, without evolution (m = 0), gen- 
eration index 7 9 = 2.7, and with a maximum generation en- 
ergy -E m ax = oo. Also shown are the HiRes observed events. 



E. Ultra High Energy Cosmic Rays Spectrum 

To accomplish the computation of the theoretical spec- 
trum we need information about the dynamical processes 
taking place in the propagation of protons along the 
CMBR. As we already emphasized, the most important 
reactions taking place in the description of a proton's 
propagation are, the pair creation p + 7 — > p + e~ +e + , 
and the photo-pion production p + 7 — > p + n. This 
last reaction is mediated by several channels. The main 
channels are 



P + l 



N + TT 
A + 7T 

R 

N + p(770) 
N + lo(782). 



(26) 
(27) 
(28) 
(29) 
(30) 



The total cross sections and inelasticities of these pro- 
cesses are well known and can be used in (|16H to compute 
the main time life of protons as a function of their energy. 
Then, with the help of expressions (|22|l and 1251) , we can 
finally find the predicted spectrum for UHECR. 

FigE] shows the obtained spectrum J{E) of UHECR 
and the HiRes observed data (two detectors, HiRes-I and 
HiRes-II). In order to emphasize the appearance of the 
GZK-cutoff in the spectrum, we have selected the ide- 
alized case when the maximum generation energy E max 
for the emitted particles from sources, is E max = 00. To 
fit the HiRes data, the generation index for the theoret- 
ical spectrum shown in the figure is 7 g = 2.7, while the 
evolution index is m — 0. Additionally, the volumetric 
luminosity is Cq = 2.96 x 10 51 ergs/Mpc 3 yr. 




10 i» 10 ia 1Q2U 1rjZ i 

Energy E [eV] 

FIG. 2: UHECR spectrum and AGASA observations. The 
figure shows the UHECR spectrum J(E) multiplied by -E 3 , 
for uniform distributed sources, without evolution, and with 
a maximum generation energy E max = 00. Also shown are 
the AGASA observed events. The best fit for the low energy 
sector (E < 4 x 10 19 eV) corresponds to 7 9 = 2.7. 



FigEl shows the obtained spectrum J(E) of UHECR 
and the AGASA observed data. Again, we have se- 
lected for the theoretical spectrum J(E) the idealized 
case -Emax = 00. To reconcile the data of the low en- 
ergy region (E < 4 x 10 19 eV), where the pair creation 
dominates the energy loss, it is necessary to have a gen- 
eration index j g — 2.7 (with the additional supposition 
that sources do not evolve) and a volumetric luminos- 
ity C — 4.7 x 10 51 ergs/Mpc 3 yr. It can be seen that 
for events with energies E > 4 x 10 19 eV, where the en- 
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ergy loss is dominated by the photo-pion production, the 
predicted spectrum does not fit the data well. To have 
a statistical sense of the discrepancy between observa- 
tion and theory, we can calculate the Poisson probability 
P of an excess in the five highest energy bins. This is 
P = 1.1 X 1CP 8 . Another statistical measure is provided 
by the Poisson x 2 given by 

x 2 = Y^i 2 ( N i h - N i hs ) + 2N ? hs H N i hs / N t)}- ( 31 ) 



Computing this quantity for the eight highest energy 
bins, we obtain x 2 = 29. These quantities show how 
far the AGASA measurements are from the theoretical 
prediction given by the curve of FigO Other more so- 
phisticated models have also been analyzed in detail 0, 
nevertheless, it has turned out that conventional physics 
does not have the capacity to reproduce the observations 
from the AGASA collaboration group in a satisfactory 
way. 

Whether HiRes or AGASA data are pointing in the 
right direction to describe the correct pattern present in 
the arrival of UHECR is still an open issue (see, for ex- 
ample, reference [43 | for a detailed comparison between 
both experimental results). In the rest of the paper we 
shall focus our attention on the possibility of an absence 
of the GZK cutoff as a consequence of LQG effects. For 
this reason, we will later return to the AGASA obser- 
vations in order to contrast the results of the following 
sections. 



III. LOOP QUANTUM GRAVITY 

Loop quantum gravity is a canonical approach to the 
problem of gravity quantization. It is based on the con- 
struction of a spin network basis, labelled by graphs 
embedded in a three dimensional insertion £ in space- 
time. A consequence of this approach is that the quan- 
tum structure of space-time will be of a polymer-like 
nature, highly manifested in phenomena involving the 
Planck scale l p . 

The above, and very brief outline of loop quantum 
gravity allows us to figure out how complicated a full 
treatment of a physical phenomena could be when the 
quantum nature of gravity is considered, even if the phys- 
ical system is characterized by a flat geometry. It is pos- 
sible, however, to introduce a loop state which approxi- 
mates a flat 3-metric on £ at length scales greater than 
the length scale C ^> l p . For pure gravity, this state is 
referred to as the weave state \W), and the length scale 
C as the weave scale. A flat weave \W) will be charac- 
terized by C in such a way that for distances d <C C the 
quantum loop structure of space is manifest, while for 
distances d ^ C the continuous flat geometry is regained. 
With this approach, for instance, the metric operator q a b 
satisfies 



A generalization of the former idea, to include matter 
fields, is also possible. In this case, the loop state repre- 
sents a matter field ip coupled to gravity. Such a state is 
denoted by \W, and, again, is simply referred to as the 
weave. As before, it will be characterized by the weave 
scale C and the hamiltonian operators H^, are expected 
to fulfill a relation analog to l|32|) . that is, we shall be 
able to define an effective hamiltonian H^, such that 



(33) 



An approach to this task has been performed by J. Alfaro 
et al. [11113113 for 1/2 -spin fermions and the electro- 
magnetic field. In this approach the effects of the loop 
structure of space at the Planck level are treated semi- 
classically through a coarse-grained approximation |45|| . 
This method leads to the natural appearance of LIV's in 
the equations of motion derived from the effective hamil- 
tonian. The key feature here is that the effective hamilto- 
nian is constructed from expectation values of dynamical 
quantities from both the matter fields and the gravita- 
tional field. In this way, when a flat weave is considered, 
the expectation values of the gravitational part will ap- 
pear in the equations of motion for the matter fields in 
the form of coefficients with dependence in both scales, 
C and l p . When a flat geometry is considered, the expec- 
tation values can be interpreted as vacuum expectation 
values for the considered matter fields. 

A significant discussion is whether the Lorentz symme- 
try is present in the full LQG theory (as in its classical 
counterpart) or not |4g. For the present work, we shall 
assume that Lorentz symmetry is indeed present in the 
full LQG theory. This assumption, jointly with the con- 
sideration that the new corrective coefficients are vacuum 
expectation values, leads us to consider that the Lorentz 
symmetry is spontaneously broken in the effective theory 
level. 

In what follows we will briefly summarize the obtained 
equations of motion for both, 1/2-spin fermions and pho- 
tons, as well as the obtained dispersion relations. 



A. Fermions 

The LQG effective equations of motion for a 1/2-spin 
fermion field, coupled to gravity, are j3^ 



d 



iAff ■ V 



B 
2C 



9 . j - „ B 



m(C - iDa ■ V)x(z) = 0, (34) 
m(C - iDa ■ V)£(x) = 0, (35) 



(W\q ab \W) =5 ab + 0(l p /£). 



(32) 



where £(x) and x( x ) are the spinor components for the 
Dirac field ^(x) = (£(x),x(x)) and the hermitian opera- 
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tors, A and C, are given by the following expressions 



the parity violation coefficients) , and where now we have 



A = 1 + + +'Jl Ijv 2 . (36) 



and the constants C and D are given by 



C = l + « 8 ^, 

2ft V 



(37) 

(38) 
(39) 



In the above expressions the m quantities are unknown 
coefficients of order 1 which need to be determined. In 
the case that $(x) is a Majorana field, the £,(x) and x{ x ) 
spinors fulfill the reality condition: 

€{x) = -ia 2 X *(x) and X (%) = ^ TO). (40) 

With the help of this condition, equations l|34l) and l|35|) 
can be simplified to 



d_ 

dt 



%B_ 
2C 



fa) 



-m(C ~iDa -\7)a 2 e{x) =0 



(41) 



Equations (|34|l and l|35|l are invariant under charge con- 
jugation C and time inversion T, but not under parity 
conjugation P. As a consequence, the fermion equation 
of motion violates the CPT symmetry through P. The 
terms that produce the P violation are those related with 
B and D. 

Some comments need to be made at this stage. Of 
specially importance to the development of the above ef- 
fective equations of motion is that they are only valid in 
a homogenous and isotropic system. From the point of 
view of a spontaneous symmetry breakup such a system is 
unique and, therefore, a privileged reference frame). It is 
possible then to put the equations of motion (and there- 
fore the dispersion relations) in a covariant form through 
the introduction of a four-velocity vector explicitly de- 
noting the existence of a preferred system. From the 
cosmological point of view, such a privileged system does 
exist, and corresponds to the CMBR co- moving reference 
system. For that reason, we shall assume that the pre- 
ferred system denoted by the presence of LIV's is the 
same CMBR co-moving reference frame, and will use it 
as the laboratory system. 

The dispersion relation for fermions can be easily ob- 
tained through the development of the Klein-Gordon-likc 
equation. The obtained dispersion relation is 

E 2 ± = (Ap±^) 2 +m 2 (C± Dp) 2 (42) 

where the ± signs correspond to the helicity state of the 
described particle (note that these signs are produced by 



1 ( ^P \ 1 K7 ;2 2 



' C 

c = 1 + K8 ^, 
D = 

2 p - 



(43) 



For our purposes, it will be sufficient to consider the lower 
contributions in both scales, l p and £ [24| 



E± = p 2 + 2ap 2 + r)p A ± 2\p + m 2 



(44) 



where we have defined the new set of corrections a, r] and 
A depending on the scales C and l p in the following way: 



a = K a (l p /C) 2 
A = k x I p /2C 2 , 



(45) 
(46) 
(47) 



being K a , k v and K\ adimensional parameters of order 1. 



B. Photons 

For the electromagnetic sector of the theory we have 
the following set of effective equations 

A(Vx5V^ + 2^ 3 V 2 (VxB) 

-26 s l p V 2 B + A9 A C 2 (Jpj l 2 p V x (B 2 B) = 0, (48) 
1 3H 

A(Vx £) + -— + 2l 2 p 9 3 V 2 {V x E) 



-29 8 l P V 2 E = 0, (49) 



where 



A=l 



2+2T 



(50) 



To calculate a dispersion relation for photons we need 
to consider only the linear part of equations l|48|l and 

(|4*9"j) . and try solutions of the type E = EQe l{ - k ' s ~ ut "> and 

B = Soe 1 ^' 1- "*', for the electric and magnetic fields. 
In this way, the obtained dispersion relation between the 
energy u> and the momentum k of photons is 



lu± = k[A 7 - 9 3 (l p k) 2 ± s l p k} 



where 



A~f — 1 + K 7 ( 



2+2T 



(51) 



(52) 
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In the previous expression the k 7 and Oi coefficients are 
adimensional parameters of order 1. As before, the ± 
signs refer to the helicity state of the described photons. 
The T quantity is a free parameter that measures a pos- 
sible non-canonical scaling of the gravitational expecta- 
tion values in the semiclassical state (let us note that the 
presence of T in the fermionic sector was not considered 
in 39] ). To be consistent with the dispersion relation of 
fermions, we shall consider only possibilities T = -1/2, 
0, 1/2, 1, etc., in such a way that Ay ~ 1 + 0[(l p /£) n ], 
where n = 2 + 2T is a positive natural number. With this 
supposition, we can find a tentative value for T, through 
the bound of the lower order correction SA ~ 0[(l p /£) n ] 
(where SA = A y — A a , being a another particle). 

Considering the lower order contributions in both 
scales, l p and £, we are able to simplify the photon dis- 
persion relation to 



k 2 + 2a 7 fc 2 ± 20 7 l p k 3 , 



where a 7 is defined by the relation A~ ( = 1 



(53) 



^7 V'Pl £) 



2+2T 



IV. THRESHOLD CONDITIONS 

A useful discussion around the effects that LIV's 
can have on the propagation of UHECR can be raised 
through the study of the threshold conditions for the re- 
actions to take place To simplify our subsequent 
discussions, let us use the following notation for the mod- 
ified dispersion relations 



E z 



P 2 + f{p) + m 2 



(56) 



where f(p) is the deformation function of the momentum 
P- 

A decay reaction is kinematically allowed when, for 
a given value of the total momentum pb = ^initial P = 
SfinaiP' one can nn d a total energy value Eq such that 
Eq ^ E m i n . Here E m i n is the minimum value attainable 
by the total energy of the decaying products for a given 
total momentum po- To find E mini it is enough to take 
the individual decay product momenta to be collinear 
with respect to the total momentum po and with the 
same direction. To see this, we can variate Eq with the 
appropriate restrictions 



Other Particles 



£ = ^£i(ft)+0(W -$>/), 



(57) 



We have so far examined the dispersion relations com- 
ing from LQG for both 1/2-spin fermions and photons. 
A relevant issue for the following developments is the 
establishment of a valid extension of the former results 
for other particles. In particular, we are interested in 
counting with dispersion relations for 3/2-spin fermions 
and 0-spin massive bosons. A precise and rigorous pro- 
cedure would require a complete calculation of effective 
field equations of motion coming from LQG for each par- 
ticle flavor in which we are interested. For present pur- 
poses we will assert that the valid dispersion relation for 
more general fermions is simply 



E\ = p 2 + 2ap 2 + rjp 4 ± 2Xp - 



m 



(54) 



This assertion preserves the basic symmetries and as- 
sumptions that have lead to the obtention of the equa- 
tions of motion for 1/2-spin fermions. 

On the other hand, in the case of bosonic 0-spin par- 
ticles, we will assert that the valid dispersion relation 
consists of 



E 2 = p 2 + 2ap 2 + T]p 4 + : 



(55) 



This assertion is based on the fact that the symmetries 
involved in the construction of the effective hamiltonian 
for 0-spin bosons would prevent the appearance of terms 
like A (which depends on the helicity). 

To conclude, let us mention that the dispersion relation 
(|54*jl will be used for the physical description of electrons, 
protons, neutrons, and A and N baryonic resonances. 
Meanwhile, the dispersion relation l|55|l will be used for 
mesons tt, p and lj. 



where £j are Lagrange multipliers, the i index specifies 
the ith particle and the j index the jth vectorial compo- 
nent of the different quantities. Doing the variation, we 
obtain 



dEj 
dpj 



— = V i 



(58) 



That is to say, the velocities of all the final produced par- 
ticles must be equal to £. Since the dispersion relations 
that we are treating are monotonously increasing in the 
range of momenta p > A, the momenta can be taken as 
being collinear and with the same direction of the initial 
quantity po . 

In this work, we will focus on those cases in which 
two particles (say a and b) collide to subsequently de- 
cay in the aforementioned final states. For the present 
discussion, particles a and b have momenta p a and pi> 
respectively, and the total momentum of the system is 
Pq. It is easy to see from the dispersion relations that we 
are considering, that the total energy of the system will 
depend only on p a = \p a \ and pb = \pb\- Therefore, to 
obtain the threshold condition for the mentioned kind of 
process, we must find the maximum possible total energy 
B max of the initial configuration, given the knowledge of 
p a and pb- To accomplish this, let us fix p a and variate 
the incoming direction of pb = npb in 

E = E a {p -p b n) + E b (pb)+x(n 2 -l)- (59) 

Varying l|59l) with respect to n (x is a Lagrange multi- 
plier), we find 



<Pb 
2X ' 



(60) 
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In this way we obtain two extremal situations \ 
±v a p b /2, or simply 



h 1 =±-a. 

v„ 



(61) 



A simple inspection shows that for the dispersion rela- 
tions that we are considering, the maximum energy is 
given by n l = —v l a /v a , or in other words, when a frontal 
collision takes place. 

Summarizing the threshold condition for a two particle 
(a and b) collision and subsequent decay, can be expressed 
through the following requirements: 



E a + E b ^J2 E f> 



(62) 



final 



with all final particles having the same velocity (vi 
for any to final particles i and j), and 



Pa-Pb 



final 



(63) 



where the sign of the momenta J2dnni Pf ^ s given by the 
direction of the highest momentum magnitude of the ini- 
tial particles. 

Our interest in the next subsections is the study of the 
reactions involved in high energy cosmic ray phenom- 
ena through the threshold conditions. To accomplish 
this goal through simple expressions that are easy-to- 
manipulate, we shall further use, for the equal velocities 
condition, the simplification 



E b m a — E a m b , 



(64) 



valid for the study of parameters coming from the region 
f(p) <C m 2 . This simplification will allow the achieve- 
ment of bounds over the order of magnitude of the dif- 
ferent parameters involved in the modified dispersion re- 
lations, which are precisely our main concern. 

In the following subsections we will study the kine- 
matical effects of LIV's through the threshold conditions 
for the reactions involved in the propagation of UHECR. 
Since, in this phenomena, photons are present in the form 
of low energy particles (the soft photons of the CMBR), 
the LQG corrections in the electromagnetic sector of the 
theory can be ignored. LQG corrections to the electro- 
magnetic sector, however, have already been studied for 
other high energy reactions such as the Mkn 501 7-rays 



the following condition must be satisfied 

m p (rn p + m 



2 5a El + l^5r] + 3i lp - 
-2E lr {\\ p \±\ p ) + AE 7T uj > 



in 

,2 



El 



m~ (2m p + m w ) 



(65) 



where E v is the energy of the emergent pion, 5 a = a p — 
0% and Srj = rj p — ry^. In expression (|65|l . the ± signs 
refer to the helicity of the incident proton. Since there 
will necessarily be a proton helicity that can minimize 
the term associated with A p and, therefore, minimize the 
energy configuration for the threshold condition, we must 
insert, in 165|) . the following equality 



2E„(\\ P \±\ P ) = 0. 



(66) 



In addition, we are assuming that the difference between 
k parameters from different particles are of order 1 (5k ~ 
1). Therefore, if not null, we can take r\ p to dominate over 
5r) in l|65l) . With these considerations in mind, we are left 
with 



2 5a El + 168r lp E* + AE^uj > 



m 2 ! (2m p + m n ) 



(67) 



Note that in the absence of LQG corrections, the thresh- 
old condition is simply 



4£Lw > 



m\{2m p + m^) 



(68) 



B. Resonant Production 7 + p — > A 

The main channel involved in the photo-pion produc- 
tion is the resonant production of the A(1232). It can 
be shown that the threshold condition for the resonant 
A(1232) decay reaction to occur, is 



25aE 2 +5 V E^ + 2 [(±) p \ p + |A A | ] E p 



AuiEp ^ m 



A 



(69) 



where E p is the incident proton energy, 5a = a p — a a and 
5t] = Tj p — t]a. Additionally, (±) p refers to the incident 
proton helicity. In the absence of LQG corrections, the 
conventional threshold condition is naturally reobtained: 



E p > 



4oj 



(70) 



A. Photo-Pion Production 7 + p — > p + tt 

Let us begin with the photo-pion production 7 + p — > 
p + 7T. Considering the corrections provided in the dis- 
persion relations (|44|l and (|55|l for fermions and bosons, 
we note that, for the photo-pion production to proceed, 



C. Pair Creation 7+p^p + e + + e 

Pair creation, 7 + p — > p + e + + e~, is greatly abun- 
dant in the sector previous to the GZK limit. When the 
dispersion relations for fermions are considered for both 
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protons and electrons, it is possible to find 



5 a 



m p + 2m e 



-E z 



. 3 m p (m p + 2m e )\ ( m. 



m p + 2m e 



E 4 



+Eu + \X e \E + -{\X p \ ± X P )E ^ m e {m p + m e ), (71) 

with Sa — a p — a e and St] = r) p — r] e . 

As in the case of photo-pion production, there will al- 
ways be an incident proton helicity which can minimize 
the inequality JHJ. Therefore, to study the production 
of the electron-positron pair under its threshold condi- 
tion, we shall set |A P | ± X p — 0. On the other hand, since 
our intention is to estimate an order of magnitude for 
the value of the diverse parameters present in the the- 
ory, let us ignore the 5rj term, since the presence of r\ p is 
of greater relevance (recall that we are considering that 
0(rf) = 0(5ri)). With these considerations, we obtain 

m R 2 3 m e 4 



5a — E' + ^-tj— E* + \X e \E 
m p 2 m p 

+Elo ^ m e (m p + m e ), 



(72) 



where we have also used m p + m e ~ m p , to simplify the 
above expression. 

Finally, if no corrections are present at all, the thresh- 
old condition would be reduced to the conventional one, 



As a last possibility, we shall also examine the bounds 
arising for the case in which no GZK anomaly (and there- 
fore no violations to the threshold and kinematics) really 
exists. Since the HiRes data have reached the ~ 1.8 x 10 20 
eV, we will consider the scenario in which no violation 
at all is confirmed by the data up to a reference energy 
E iei = 2 x 10 20 eV. 

In order to study the different corrections, given that 
we don't have a detailed knowledge of the deviation pa- 
rameters, we shall take account of them independently. 
Naturally, there will always exist the possibility of hav- 
ing an adequate combination of these parameter values 
that could affect the threshold conditions simultaneously. 
However, as shall soon be evident, each one of these pa- 
rameters will be significant at different energy ranges. 



1. a Correction 

We shall begin our analysis with correction a and the 
consideration of the threshold condition for pair produc- 
tion. In this case we have 



5a — -E + Eu> m e (m v + m e ), 



(74) 



with 5a = a p — a e . As is clear from the above condi- 



tion, the minimum soft-photon energy u> n 
production to occur, is 



for the pair 



Elo m e {m p + m e ). 



(73) 



D. Bounds 

In order to study the threshold conditions 1)67(1 . Iti9|) 
and H72J1. in the context of the GZK anomaly, we must 
establish some criteria. 

Firstly, as we have seen in section [HI the convention- 
ally obtained theoretical spectrum provides a very good 
description of the phenomena up to an energy ~ 4 x 10 19 
eV. The main reaction taking place in this well described 
region is the pair creation j+p — > p + e + + e~ and, there- 
fore, no modifications are present for this reaction up to 
~ 4 x 10 19 eV. As a consequence, and since threshold 
conditions offer a measure of how modified kinematics is, 
we will require that the threshold condition (|72|1 for pair 
creation not be substantially altered by the new correc- 
tive terms. 

Secondly, we have the GZK anomaly itself, which we 
are committed to explain. Since for energies greater than 
~ 8x 10 19 eV the conventional theoretical spectrum does 
not fit the experimental data well, we shall require that 
LQG corrections be able to offer a violation of the GZK- 
cutoff. The dominant reaction in the violated E > 8 x 
10 19 region is the photo-pion production and, therefore, 
we shall require further that the new corrective terms 
present in the kinematical calculations be able to shift 
the threshold significantly to preclude the reaction. 



m e m e 
w min = -^r(m p + m e ) - da — E. 
E m p 



(75) 



It follows therefore that the condition for a significant 
increase or decrease in the threshold energy for pair pro- 
duction becomes |<5a| ^ m p (m p + m e ) / E 2 . In this way, if 
we do not want kinematics to be modified up to a refer- 
ence energy E le { = 3 x 10 19 , we must impose the following 
constraint 



(m p + m e )m p 
E± t 



= 9.8 x 10" 



(76) 



Similar treatments can be found for the analysis of other 
astrophysical signals like the Mkn 501 7-rays when 
the absence of anomalies is considered. 

Let us now consider the threshold condition for the 
photo-pion production. Taking only the a correction, we 
have 



2 5a El + AE^uj ^ 



to 2 (2to p + m n ) 
m p + m-a 



(77) 



It is possible to find that for the above condition to be 
violated for all energies E w of the emerging pion, and 
therefore no reaction to take place, the following inequal- 
ity must hold 
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where uj = KT = 2.35 x 10~ 4 eV is the thermal CMBR 
energy. If we repeat these steps for the A(1232) resonant 
decay, we obtain the following condition 

a A -a p > — = 1.7 x 1(T 25 [lj/oj ] 2 . (79) 

m\ — m£ 

To estimate a range for the weave scale C, let us use 
as a reference energy w rc f = w m i n , where w m i n is the min- 
imum energy for the reaction to take place, in inequal- 
ity 177(1 . when the condition for a significant increase in 
the threshold condition is taken into account (for a pri- 
mordial proton reference energy E rc { = 2 x 10 20 , this is 
^min ~ 2.9 x wo) 7 and join the results deduced from the 
mentioned requirements. Assuming that the n a param- 
eters are of order 1, as well as the difference between 
them for different particles, we can estimate — for the 
weave scale C — the preferred range 

2.6 x 1(T 18 cV" 1 < C < 1.6 x 1CT 17 eV -1 , (80) 

where the lefthand and righthand sides come from 
bounds (|76|l and (|78|l respectively (since the A(1232) is 
just one channel of the photo-pion production, we shall 
not consider it to set any bound). 

If no GZK anomaly is confirmed in future experimental 
observations, then we should state a stronger bound for 
the difference a n — a p . Using the same assumptions to set 
the restriction (|76l) when the primordial proton reference 
energy is E le { = 2 x 10 20 eV, it is possible to find 

K - a p \ < 2.3 x 1CT 23 . (81) 

In terms of the length scale C, this last bound may be 
read as 

C > 1.7 x KT 17 eV~\ (82) 

which is a stronger bound over C than l|76|) . offered by 
pair creation. 



which is, of course, a strong bound over a parameter of 
order 1. 

Since the basis of the effective LQG methods which we 
have developed rely on the fact that the coefficients k are 
of order 1, we must conclude that a reaction of the type r\ 
should be discarded, in opposition to the expectations of 
our previous work [24j , when only photo-pion production 
was analyzed. 



3. A Correction 

Finally we have to consider the A correction. In our 
previous work |24j . having studied the photo-pion pro- 
duction through the A(1232) channel decay, we empha- 
sized the possibility of a helicity dependent violation. For 
this effect to take place, the following configuration must 
be satisfied (using again u> rc { = w m i n ) 

|A P | > |A A | + 1.3 x 1CT 3 eV. (86) 

However, when pair production is analyzed (in the same 
way as with a and ry), the following condition emerges 
for A e : 

|A e | < 1.6 x 1(T 5 eV, (87) 

which is more than 1 order of magnitude stronger than 
the required value for protons in (|86|l . This weakens the 
possibility of a limit violation through helicity dependent 
effects. 

Another stronger bound can be found in psj , where a 
dispersion relation of the type 

E 2 = p 2 + Xp + m 2 (88) 

is analyzed, and it is found that the case A ^ 1 x 10~ 7 
eV should be discarded because of the highly sensitive 
measurements of the Lamb shift. 



2. n Correction 

Let us now turn our attention to the rj parameter. The 
threshold condition for the pair production, when only 
the rj parameter is considered, is 

—rj — -E 4 + Euj ^ m e (m p + m e ). (83) 

Repeating the same analysis we did for the a parameter, 
it is possible to find the following constraint 

M < ;^rK + m e ) = 1.6 x 1CT 60 eV- 2 . (84) 

Recalling that rj = n v l p , result l|84|l can be reexpressed 
in the form 

\n n \ < 2.4 x KT 4 , (85) 



E. The Cubic Correction 

A commonly studied correction which has appeared in 
several recent works |23l H3| , and which deserves our 
attention, is the case of a cubic correction of the form 

E 2 =p 2 +m 2 +£p 3 , (89) 

(where £ is an arbitrary scale). It is interesting to 
note that strong bounds can be placed over deformation 
f(p) = £p 3 . We will assume in this section, that £ is an 
universal parameter (an assumption followed by most of 
the works in this field). 

^From dispersion relation l|89|) , the threshold condition 
for the photo-pion production is 

4w£ + 2£- p , n E 3 ^ mJ2m„ + mJ. (90) 

(m p + m^Y 
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Clearly, when £ is negative, it can be observed from l|90|) 
that the threshold energy for photo-pion production can 
be easily shifted, preventing the reaction to take place 
for high energies. The condition for this to be the case is 



128 



27 m 3 m p \2m p + m 
= 7.75 x 1CT 45 [w/w ] 3 eV~ 



(91) 



It can be seen therefore that in the particular case of 
£ = — l p (lp = 8.3 x 10~ 29 eV _1 ), the reaction can be 
considerably suppressed. 

Let us note, however, that the pair production p + 7 — > 
p + e + + e~ imposes strong restrictions over a negative £ 
parameter. Following the methods that we have used up 
to now, it is possible to find that the threshold condition 
for the pair production is 



(m p + 2m e ) 2 



-E > m P 



(92) 



Since we cannot infer modifications in the description of 
pair production in the cosmic ray spectra up to energies 
£~4x 10 19 eV, we must at least impute the following 
inequality 



\i\ < 



(m p + m e )(m p + 2m e Y 
m p Ef c{ 

= 3.26 x 1CT 41 eV" 1 
(= 3.98 x HT 13 ^), 



(93) 



where we have used E IC { = 3 x 10 eV. This last result 
shows the strong suppression over £ . As a consequence, 
the particular case |£| = l p = 8 x 10~ 28 cV _1 should be 
discarded. A bound like (|93|l seems to have been omitted 
up to now for the GZK anomaly analysis. 



V. MODIFIED SPECTRUM 

In this section we shall show how the only surviving 
LQG correction from our previous analysis in section llVl 
A = 1 + a, can affect the prediction of the theoretical 
cosmic ray spectrum. Our approach will be centered on 
the supposition that the LQG corrections to the main 
quantities for the calculation — such as cross sections and 
inelasticities of processes — are, in a first instance, kine- 
matical corrections, and that the Lorentz symmetry is 
spontaneously broken. These assumptions will allow us 
to introduce the adequate corrections when a modified 
dispersion relation is known. 



A. Kinematics 

When having a spontaneous Lorentz symmetry break- 
ing, we can use the still valid Lorentz transformations to 



express physical quantities observed in one reference sys- 
tem, in another one. This is possible since, under a spon- 
taneous symmetry breaking, the group representations of 
the broken group preserves its transformation properties. 
In particular, it will be possible to relate the observed 4- 
momenta in different reference systems through the usual 
rule 



P» = A »P»> 



(94) 



where p^ = {—E,p) is an arbitrary 4- momentum ex- 
pressed in a given reference system /C, p'^ is the same vec- 
tor expressed in another given system fC' , and A " is the 
usual Lorentz transformation connecting both systems. 
Such a transformation will keep invariant the scalar prod- 
uct 



-E 2 



P 



(95) 



as well as any other product. 

Let us illustrate, for transformation Q94|). the situation 
in which KJ is a reference system with the same orienta- 
tion of K, and which represents an observer with velocity 
(3 with respect to K.. In this case shall correspond to 

a boost in the (3 = (3/\(3\ direction, and expression 194|) 
will be reduced to 



E' = j(E-p-p), 
p' = 1{P-$E), 



(96) 
(97) 



where 7 = (1 — 1 ) 1//2 . A particular case of this trans- 
formation will be that in which has the same direction 
as p, and /C' corresponds to the cm. reference system, 
that is to say, the system in which p' — O. In such a case 
we will have (3 = p/E and 7 = E/ (E 2 — p 2 ) -1 / 2 , jointly 
with the relation 



E' = E/j = (E 2 -p 2 ) 



2\1/ 2 



(98) 



In other words, the cm. energy of a particle with energy 
E and momentum p in /C will correspond to the invariant 
(E 2 — p 2 ) 1 / 2 . Furthermore, such energy is the minimum 
measurable energy by an arbitrary observer; this can be 
confirmed by solving equation dE' /d(3 = from the re- 
lation (|96|l and by verifying that the solution is (3 = p/E. 
This allows us to interpret E' — (E 2 — p 2 ) 1 / 2 as the rest 
energy of the given particle. To simplify the notation 
and the ensuing discussions, let us introduce the variable 
s = (E 1 ) 2 , where E' is given by (|9"5|l . 

So far in our analysis, relativistic kinematics has not 
been modified. Nevertheless, a difference with the con- 
ventional kinematical frame is that in the present theory 
the product l|95|l will not be independent of the particle's 
energy; conversely, we will have the general expression 



P^Pfj, 



-fa(E,p) 



(99) 



where f a (E,p) is a function of the energy and the mo- 
mentum, that represents the LIV provided by the LQG 
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effective theories. Let us note that expression (|99|l is just 
the modified dispersion relation 



E 2 



-p 2 + f a (E,p)+m 2 a 



(100) 



To be consistent f a (E,p) must be invariant under 
Lorentz transformations and, therefore, can be written 
as a scalar function of the energy and the momentum. 

We have already made mention of the fact that LIV's 
inevitably introduce the appearance of a privileged sys- 
tem; in the present discussion we will choose as such a 
system the isotropic system (which by assumption is the 
co-moving CMBR system), and will express f a (E,p) in 
terms of E and p measured in that system. As may be 
expected in this situation, f a will be a function only of 
the energy E and the momentum norm p = \p\ , since no 
trace of a vectorial field could be allowed when isotropy 
is imposed. For example, in the particular case of the 
dispersion relation for a fermion, the function f a (E,p) 
depends uniquely on the momentum, and can be written 
as 



f a {p) = 2a a p 2 + r] a p 4 ± 2X a p. 



(101) 



For simplicity, we shall continue using f a {p) instead of 

fa(E,p). 

Through the recently introduced notation and the use 
of expression (|98|l . the cm. energy of an a particle with 
mass m a and deformation f a (p) will be 



# 3 = Vfa(p)+m 2 . 



(102) 



Of course, the validity of this interpretation will be sub- 
ordinate to those cases in which 



fa(p) + m 2 a > 0, 



(103) 



or, equivalent, to those states with a time-like 4- 
momentum. In the converse, particles with energies and 
corrections such that s a — fa(p) + rn 2 ^ will be de- 
scribed by light-like physical states if the equality holds, 
or space-like physical states if the inequality holds. 

A new effect provided by LIV's is that, if a reference 
system where p = exists, then in that system the par- 
ticle will not be generally at rest. To understand this it 
is sufficient to verify that in general the velocity follows 



BE p 
~3p^E' 



(104) 



and therefore does not vanish at p = 0. Returning to 
equations (|96|) and 1)9 7|). we can see that when (3 — v — 
BE /dp, the following result is produced 



BE' 
Bp' 



7« 



BE 
Bp 



Bp 
Bp' 



0. 



(105) 



where j v = (1 — v 2 )^ 1 / 2 . Result 1|105|) shows that the ve- 
locity of the system where the particle is at rest is effec- 
tively v. There emerges, then, an important distinction 



between the phase velocity p = p/E of the cm. system 
of a particle, and the group velocity v = BE /Bp of the 
same particle. 

The above results, for a single particle, can be eas- 
ily generalized to a system of many particles. For in- 
stance, the total 4-momentum of a system of many par- 
ticles p^ Qt = Pi wm transform through the rule 194|l , 
and the scalar product (p^Pfi) tot will be an invariant un- 
der Lorentz transformations (as well as any other prod- 
uct). As in the case of individual particles, we define y/s 
as the total rest energy measured in the system of cm. 
That is to say: 

s = E 2 ot -p 2 ot . (106) 

In the case in which we have a system composed of a 
proton with energy E and momentum p, and a photon 
(from the CMBR) with energy u> and momentum k (all 
these quantities are measured in the laboratory isotropic 
system /C), the s quantity shall acquire the form 

s = {E + oj) 2 - (p + k) 2 

= {E' + u'f- (p' + k) 2 , (107) 

where the E' and p' quantities are measured in an ar- 
bitrary reference system. In particular we are interested 
in the system where the proton momentum p' is null; 
that is to say, the system in which E' — ^ fsZ . If e is the 
photon energy in such a system, then 



= + (108) 

where we have used the dispersion relation ui — k (or 
e = k') for the CMBR photons. These results will allow 
us to express the main kinematical quantities in term of 
e and s p = f p (p) + m 2 . For example, the Lorentz factor 
that connects the JC system with the JC' system where 
p' = 0, will be 



E 



7 



(109) 



Meanwhile, the Lorentz factor connecting JC with the cm. 
system (that in which p' + k = 0), will be 

E + uj 

7^ = T 

s P + 

~ § . (HO) 

As a last comment, let us note that to the first order in 
the expansion of the dispersion relations in terms of the 
scales C and l p , when we consider high energy processes 
such that p 2 ^> f(p) + m 2 we can freely interchange the 
momentum p by the energy E in the deviation function 
f(p). That is to say, we may consider as a valid relation 
the following expression 

E 2 =p 2 + f(E)+p 2 , (111) 

where we have made the replacement f(p) — > f(E). This 
procedure will greatly simplify the next discussions. 
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B. Modified Inelasticity: p + 7 — > p + x 

Following the same methods of section^J let us obtain 
a modified inelasticity K for a process of the type p+7 — > 
p + x, where x is an emitted particle that, in the present 
physical problem in which we are interested, can be a w, 
p or lu meson. We note that the dispersion relation for 
the emerging proton (after a collision with a photon) can 
be written in the form: 

E 2 ~p 2 p = f p (E p )+m 2 p , (112) 

where E p is the final proton energy. Since the left side of 
(|112f> is invariant under Lorentz transformations, we can 
write 

(E;) 2 -(p* p ) 2 = f p (E p ) + m 2 pl (113) 

where the * denotes the quantities measured in the cm. 
system. On the other hand, in such a system, the follow- 
ing conservation relations of energy and momentum are 
satisfied: 

e; + e* x = v~s (in) 

and 

(p;r = (pd 2 - (us) 

Substituting both quantities in relation l|112l) . we can 
obtain 

2^E* p = s + f p (E p ) - f x (E x ) + m\ - m 2 x , (116) 

or, in a more convenient form 

2y/EE; = s + s p (E p ) - s x (E x ). (117) 

In the same way, we also have the energy conservation 
relation in the laboratory system: 

E p + E x = E tat . (118) 

Using the definition for the inelasticity K x = AE/E for 
a process, where AE = Ei — Ef ~ E tot ~Ef, it is possible 
to rewrite l|118|) in terms of K x through expressions 

E x = K X E, (119) 
E p = (1-K X )E, (120) 

where E is the initial energy of the initial proton. Having 
done this, equation (|117H now acquires the form: 

2yfsE* p = s + s p [(l - K X )E] - s x [K x E\. (121) 

To simplify the development of the inelasticity, let us 
write the former relation as E* = F(E,K X ), where F — 
F(E,K X ) is defined through 

F = -L(s + s p [(l - K X )E] - s x [K x E}). (122) 

l\l s 



On the other side, the Lorentz transformation rules give 
us the relation between the proton energies in the lab- 
oratory system and the cm. system respectively. This 
relation is 

E p = lc (E* p +f3 cP * p cos 9) 

= 7c (E* p + p c y/E;* - s p (E p ) cosfl) . (123) 

Joining (|122fl and l|123|) . it is possible to find the general 
equation for K x : 

{l-K x )yTs=(F(E,K x ) 

+ ^F 2 (£, K x ) - Sp [(l - K X )E] cos^) . (124) 

It should be noted, however, that the solution for K x 
from l|124(l will depend in the 9 angle. For this reason, 
once this last equation is resolved, it is convenient to 
define the total inelasticity K as the average of K x with 
respect to the 9 angle. That is to say 

1 r 

K = - K x d9. (125) 
n Jo 

It is relevant to mention that now, as opposed to result 
IjlUp. the inelasticity K will be a function of both, the 
energy E of the initial proton and the energy e of the 
CMBR photon. 

C. The m\ — > s a = m 2 a + f a (E) Prescription 

Let us recall our interpretation relative to the fact that 

Sa" 2 = (fa{E a ) + m 2 ) 1 ^ 2 can be understood as the rest 
energy of a particle a, as a function of the energy E a 
that it has in the laboratory system JC. As we have al- 
ready emphasized, such an interpretation will be valid 
for particles with time like 4-momenta. 

In the reactions given between high energy protons and 
the photons of the CMBR, the whole scenario consists of 
the collision between two particles, p and 7, with the sub- 
sequent production of a certain number of final particles. 
Let us suppose that a is one of these particles in the final 
state. The knowledge of the inelasticity K for the reac- 
tion will allow us to estimate the average energy < E a > 
with which such a particle emerges (since K provides the 
average fraction of energy with which such a particle is 
produced). That is to say, on average, the rest energy of 
the final particle a will be s]J 2 — [f a {< E a >) + m 2 } 1 / 2 . 
Moreover, the knowledge of the inelasticity K will allow 
us to express s a as a function of the energy E of the 
initial proton: 

s a = s a (E). (126) 

Following our previous interpretation, we can view the 

recently described process as a reaction between a pro- 

1/2 

ton with mass s p , which loses energy emitting particles 
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a with mass s~J 2 calculated in the previous form. This 
idealized reasoning gives us a clear prescription to kine- 
matically modify those dynamical quantities with which 
we must work and where energy conservation is involved. 
This prescription is: 



ml->s a (E) = f a (E)+: 



(127) 



where we have expressed correction /„ as a function of 
the initial energy of the incident proton. 

Prescription Ijl27(l establishes the notion of an effective 
mass which is dependent on the initial energetic content 
of a reaction. As a consequence, given the explicit knowl- 
edge of the dependence that a cross section has on the 
masses and energies of the involved states, to obtain the 
modified version, it will be appropriate to use the dis- 
cussed prescription. 

Let us note, however, that an important weakness of 
the present method is the inability to determine whether 
nip comes from the initial state proton (with an initial 
energy E), or the final state proton (with a final energy 
E p = (1 — K)E); specially since that distinction gives 
rise to different values for s p . To overcome this difficulty, 
and therefore any ambiguity in the prescription, we shall 
restrict our treatment to the case a„ = 0. 



D. Redshift 

Another important problem related to the introduction 
of LIV's in the dispersion relations is whether the red- 
shift relation for the propagation of particles in a FRW 
universe is modified. This could be of great relevance be- 
cause of the large distances involved in cosmic ray propa- 
gation and, therefore, the possible cumulative effects. We 
shall examine this issue through the study of a classical 
point-like particle propagating in a FRW space-time. 

The most general action for a point particle in a given 
space-time is 



J a 



Adr, 



(128) 



with t an affine parameter chosen to accomplish dr = 
J— g^udx^dx" (we are using the mostly plus signature 
(— , +, +, +)). The variation of the action can be realized 
through two separated terms: 



SS 



b 

ASdr + / 5 A dr. 



(129) 



Starting with the first term, it is possible to write Sdr 
in the following way (using — dx^ /dr and identifying 
the free torsion connections T"); 

dr^ 

Sdr = -g^—5dx v -T v fl(7 u l *u u dT8x' T . (130) 

Naturally, in the above expression we can allow Sdx 11 = 
d(Sx fJ ') and Sdr — d(Sr). The variation of the second 



term in i|129|) can be developed through 



SA 



dA 
~dx° 



5x° 



dA 
dw 7 



Su° 



(131) 



For the <5u M variation we must proceed carefully since 
there will be constraints between 5m m and Sx^. We have 
5u^ = u'^ - u", where u'" = dx'^/dr', with dx^ = 
dx^ + Sx^ and dr' = dr + St. In this way, to the first 
order in the variations, it is possible to find the following 
constraint between Su^ and Sx^ (where we have used 
equation 113U|) to work out the relation): 



(Su^dr = d(6x' i )+ui'\g vv u r >5dx u 
+Ff la u v drSx a ]. 



(132) 



Using the former relations, the complete expression for 
the variation of the action is 



5S = 



OA 



Sx" 



b r 



dA 

dx^ 



. _ — yiy* ,, d ( dA 



where we have defined A = (OA/ du a )u° 
define the momenta p^ as follows 



SS 
8x^ 



dA 



drSx^, (133) 
- A. Let us now 

(134) 



This definition is concomitant with the canonical ap- 
proach and allow us to write the equations of motion 
in a very simple and convenient way: 



dPt, = 



dA 

dx^ 



Ar>^ 



dr. 



(135) 



We are interested in obtaining an expression for the 
redshift relation having as a starting point the equations 
of motion H135fl deduced from the S action. For this we 
must consider the FRW metric 



g^ = diag 



-1, 



R 2 (t) 

1 — kr 2 



,R 2 (t)r 2 ,R 2 (t)r 2 sin 2 9 



.(136) 



To accomplish our goal, let us calculate the variation of 
P 2 = g^PiPj through a path parameterized by r: 

= (-^9 ij )PiPj + 2g ij P i^r. (137) 

Note that inn the preceding expression, dpj is given by 
the dynamical equations (|135fl . Using these equations to 
simplify equation l|137|) . it is possible to deduce that 



dr 



(pR) = -g i lp i n 



(138) 
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where fL- is defined through 



dA 

dx l 



v du k 



(139) 



and the T k , are the FRW spatial connections given by 



9 ' (9il,j +9jl,i ~ 9ij,i)- 



(140) 



To conclude, we now turn our attention to fij. If we 
want to be loyal to the spirit of the FRW space-time 
formulation, we must impose isotropy and homogeneity 
on the Lagrangian A, when expressed in the co-moving 
FRW frame (a characteristic present in our previous de- 
velopment of LIV's). That is to say, A = A(u 2 , t), where 



9ijt 



In this way, the spatial dependence will be 



through Qij . If we differentiate A with respect to x then 



dA 

dx k 



dA du 2 

du 2 dx k 
OA • j dg l] 



du 



-u u J 



dx k 



which implies that 



OA 

dx k 



dA 

du 1 



T'- 



3k- 



(141) 



(142) 



Note that this in turn will mean that fii = 0. So, as a 
general result, the usual redshift relation is reobtained: 



— (pR) = 0. (143) 



E. Spectrum and Results 

Introducing the above modifications to the different 
quantities involved in the propagation of protons (like 
the cross section a and inelasticity K), we are able to 
find a modified version for the UHECR energy loss due 
to collisions. Since the only relevant correction for the 
GZK anomaly is a, we focused our analysis on the par- 
ticular case f(p) = 2ap 2 . To simplify our model we re- 
stricted our treatment to the case a > (consistent with 
the effective mass interpretation) and used only a m =/= 0, 
where a m is assumed to have the same value for mesons 
7r, p and ui. 

FigElshows the modified energy loss t(E) for UHECR 
obtained for different values of a m . These are, curve 
1: a m = 9 x 1CT 23 (£ ~ 8.6 x 1CT 18 eV" 1 ); curve 2: 
a m = 5 x 1CT 23 (£ ~ 1.2 x KT 17 eV" 1 ); and curve 
3: a m — which corresponds to the case without mod- 
ifications given by the conventional theory. It can be 
seen therefore how the corrections can affect the main life 
time of protons propagating through the CMBR, allow- 
ing a strong improvement in the distances that protons 
can reach before loosing their characteristic energy (for 
energies greater than 1 x 10 20 eV). The effects that the 
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10 13 10* 
Energy E [eV] 

FIG. 3: Modified energy loss for UHECR due to collisions. 
The figure shows the case a m 7^ 0, for three different values of 
the weave scale £ Curve 1: a m = 9 X 10~ 23 (£ ~ 8.6 X 10" " 
eV^ 1 ); curve 2: a m = 5 x 1(U 23 (£ ~ 1.2 x 10" 
curve 3: a m = (without modifications). 



eV- 1 ); 



LQG corrections have on the propagation of UHECR are 
manifest through a decay of the energy loss in the range 
£~ lx 10 20 eV. To understand this, recall relation 1(77)1 
for the threshold condition of photo-pion production: 



2 5a El + AE^uj ^ 



to 2 (2m p + m n ) 



(144) 



As we saw in section llVl the condition for a significant 
increase or decrease in the energy threshold can be cal- 
culated as \Sa\ > {2m p + m- 7l ){mp + m- n )/2E 2 . Therefore, 
for a given value of 5a > 0, the energy at which the LIV 
effects start to take place is 



E = — — (2m p + m 7 r)(m p + m OT ). 
Ida 



(145) 



In the case a m = 9 x 10~ 23 (curve 1 of FigOJ, this energy 
is E = 1.1 x 10 20 eV, while in the case a m = 5 x 1CU 23 
(curve 2) this corresponds to E = 1.5 x 10 20 eV. Beyond 
these energy scales, at about E ~ 2 x 10 20 eV, a sharp de- 
cay is observed in the behavior of the curve. This is due 
to the fact that the modified inelasticity K will strongly 
constraint the energy-momentum phase space accessible 
to the final states depending on the initial energy E that 
the primary proton carries (recall that now K is a func- 
tion of the energy E of the incident proton and the energy 
e of the CMBR photon). 

We can also find the modified version of the UHECR 
spectrum for a m ^ 0. Fig0] shows the AGASA obser- 
vations and the predicted UHECR spectrum in the case 
a m = 1.5 x 1CT 22 (£ ~ 6.7 x 1(T 18 eV" 1 ) for three differ- 
ent maximum generation energies E max . These are, curve 
1: 5 x 10 20 eV; curve 2: 1 x 10 21 eV; and curve 3: 3 x 10 21 
eV. The Poisson probabilities of an excess in the five high- 
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FIG. 4: Modified UHECR spectrum and AGASA observa- 
tions. The figure shows the modified spectrum J(E) mul- 
tiplied by E 3 , for uniform distributed sources and without 
evolution, for the case a m = 1.5 x 1CT 22 (£ ~ 6.7 x 10~ 18 
eV" 1 ). Three different maximum generation energies E max 
are shown. These are, curve 1: 5 x 10 20 eV; curve 2: 1 x 10 21 
eV; and curve 3: 3 x 10 21 eV. 



est energy bins for the three curves are P\ = 3.6 x 10 , 
P 2 = 2.6 x 1CT 4 and P 3 = 2.3 x 1CT 4 . The Poisson X 2 
for the eight highest energy bins are X\ — 10, X2 = 10-9 
and xi = 11.2 respectively. The possibility of reconcil- 
ing the data with finite maximum generation energies is 
significant given that conventional models require infinite 
maximum generation energies E max for the best fit. For 
the lower part of the spectrum (under E = 4 x 10 19 eV), 
the parameters under consideration leave the spectrum 
completely unaffected. This is due to the fact that in 
such a region the dominant reaction is the pair produc- 
tion, which has not being modified to obtain the spec- 
trum. A more accurate study on this issue would require 
the computation of a modified inelasticity for the pair 
creation. Meanwhile, we must content ourself with the 
semiqualitative criteria given in section llVl to rule out the 
parameters. 



VI. CONCLUSIONS 

The scientific challenge that represents the search for 
new empirical backgrounds to test quantum gravity the- 
ories is at the embryonic stage. In this context, the pos- 
sibility that ultra high energy cosmic rays could be ex- 
periencing quantum gravity effects places us in a very 
challenging situation which deserve attention. Neverthe- 
less, the present stage of UHECR observations demands 
that we proceed with caution and patience. 

We have seen how the kinematical analysis of the differ- 
ent reaction taking place in the propagation of ultra high 
energy protons can set strong bounds on the parameters 



to the theory. In comparison with our previous work, we 
have eliminated some previously open possibilities by the 
particular study of the pair creation p + 7 — > p + e + + e~ , 
in the energy region where this reaction dominates the 
proton's interactions with the CMBR. In this way, the 
only possibility still open (for the corrective terms con- 
sidered in the expansion (|54|l for the dispersion relations) 
and favored by the LQG scales, is the correction a. If 
this is the case, a favored region for the scale length £ 
estimated through the threshold analysis would be 

2.6 x 1(T 18 eV -1 < £ < 1.6 x 1CT 17 eV -1 . 

Similarly, the kinematical corrections can be studied 
in more detail when their effects are considered in the 
theoretical spectrum. In this regard, we have seen how 
to develop a modified version of the inelasticity for the 
photo-pion production, and its implications in the mean 
life time of a high energy proton as well as on the spec- 
trum. To accomplish this last task we have only assumed 
a spontaneous Lorentz symmetry break up in the effective 
equations of motion, allowing the use of Lorentz trans- 
formations on the dispersion relations. Therefore, result 
(112411 can be used in a more general context than the 
special case offered by the LQG framework. 

Special mention must be made of a recent develop- 
ment |4~i| , where the dispersion relations for fermions and 
bosons are generalized to include the extra factor T. It 
should be noted that this new factor always appears in 
the dispersion relations in the form (Z p /£) T , such that 
the parameter a in our equation 144f> gets an extra factor 
\l p /£) ■ This freedom can be used to move the scale 
£ down if needed, so that the cosmic ray momentum p 
always satisfies the bound p £ 1 , without changing our 
prediction of the UHECR spectrum. 

Future experimental developments like the Auger ar- 
ray, the Extreme Universe Space Observatory (EUSO) 
and Orbiting Wide-Angle Light Collectors (OWL) satel- 
lite detectors, will increase the precision and phenomeno- 
logical description of UHECR. On the more theoretical 
side, progress in the direction of a full effective theory, 
with a systematic method to compute any correction with 
a known value for each coefficient, is one of the next steps 
in the "loop" quantization programme 49, 50]. There- 
fore, it is important to trace a phcnomcnological under- 
standing of the possible effects that could arise as well as 
the constraints on LQG, in the high and low energy regi- 
mens (for other phenomenological studies of LQG effects, 
see for example [H]J and [52^1. 
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